Entanglement distillation using particle statistics 
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We extend the idea of entanglement concentration for pure states(Phys. Rev. Lett. 88, 187903) 
to the case of mixed states. The scheme works only with particle statistics and local operations, 
without the need of any other interactions. We show that the maximally entangled state can be 
distilled out when the initial state is pure, otherwise the entanglement of the final state is less than 
one. The distillation efficiency is a product of the diagonal elements of the initial state, it takes the 
maximum 50%, the same as the case for pure states. 
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Entanglement shared between two distinguishable par- 
ticles is generally an advantage for quantum information 
processing, while indistinguishability prevents us from 
addressing the particles separately that seems to be a 
disadvantage in information processing. It is only re- 
cently that researchers have started to investigate the 
use of particle statistics (both bosonic and fermionic) for 
quantum information processing^, 0, 0- It was shown 
that quantum statistics can lead to an effective interac- 
tion between internal and external degrees of freedom of 
particles lj, and then result in space-spin entanglement 
transfer. The other useful tasks such as entanglement 
concentration and state discrimination Q could also be 
accomplished, even if nonoptimally, using only the ef- 
fects of quantum statistics, without the need of any other 
interactions. These proposals purely based on particle 
statistics differ significantly from some previous sugges- 
tions, where particle statistics together with interparticle 
interactions 0, 0, @ were used for quantum information 
processing. 

The entanglement distillation/concentration is essen- 
tial in quantum communication and computation, it in- 
creases the entanglement shared between distant par- 
ties, and consequently could better the performances 
in quantum information processing. As entanglement 
can not be increased by local operations and classical 
communication [3, the distillation/concentration opera- 
tion is to distill/concentrate entanglement loc ally from a 
largeto a smaller number of pairs [IUDU 13 III 13 13 

1313 E2. 

Quantum information processing based on particle 
statistics is very useful for tasks implemented with iden- 
tical particles. In this paper, we extend the protocol 
proposed in Ref. 2] for pure states to mixed states. To 
start with, we recall and summarize the entanglement 
concentration protocol using only particle statistics as 
follows. The initial state is prepared at two sides labelled 
by L and R. In each side, two parties Alice (A) and Bob 
(B) share n pairs of particles. The protocol includes the 
following four steps. (1) Bring each of Bob's particles 
(from both sides of L and R) into a 50/50 beam split- 
ter; (2) Before the particles from the side L fall into the 
beam splitter, a specific unitary transformation that flips 
the particle's spin is applied to each of the particles; (3) 
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FIG. 1: An schematic illustration of the scheme. Alice and 
Bob start with two pairs of entangled systems at side L and 
R, respectively. At each side, Alice and Bob share n pairs of 
particles. After a set of local operations including spin flip and 
50/50 beam splitter, path measurements, and classical com- 
munications, states with more entanglement can be distilled 
out from the initial state. The scheme works with particle 
statistics solely, without the need of other interactions. 



Make a specific measurement on the output particles, and 
discard those which bunch(for fermions)/antibunch(for 
bosons); (4) Repeat the above three steps for all par- 
ticles with Bob. This was schematically illustrated in 
figure n We now apply the above steps to mixed states 
for fermions, the results can be straightforwardly gen- 
eralized to the case of bosons. As the authors did in 
Ref. [3 , we consider the entanglement in the internal de- 
grees of freedom of the particles, for example the spin in 
the case of fermions. In either sides L or R, Alice and 
Bob share n pairs of particles, the initial state at side / 
(/ = R,L), distributed between Alice and Bob, in basis 
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where |^4t)£ = I TT ••• T)l.a denotes n spin-up particles 

n 

at L side with Alice. The other notations in the basis are 
similar to |^4f)2- The total initial state under considera- 
tion is then 
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Here, the basis chosen is 
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With this total initial state, next step we let the 
first pair of particles (labelled by 1 in the circle) with 
Bob go to the beam splitter 1, after having the spe- 
cific unitary transformation in step (2), the total state 
still takes the form of Eq.(0), but was written in ba- 
sis of {\A])l\A]) R \B])l\B\) R , \A])l\A[)l\B])l\B]) R , 
\A\)l\A]) R \B\)l\Bl)% \Al)l\A\) R \Bi)l\B]) R .} Sim- 
ilarly, the state of particles after the first pair having 
passed through the 50/50 beam splitter can be written in 
the same form as in Eq.@ in the Hilbcrt space spanned 
by 
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where \Bx) LlR1 = \i{\B ]) L1 \B i) L1 + \B ]) RX \B + 
\(\B ]) Ll \B |) R1 + \B t)m), and \B 2 ) L1R1 = 

-\i{\B T)li|S + \B T)m|5 [) R1 ) + [) Ll \B ] 
)m + \B T)li|-B 1)ri)- When the particle had passed 
through the beam splitter, Bob performs a path mea- 
surement on the first pair with assumption that the de- 
tectors do not absorb the particles and do not disturb 
their internal(spin) degrees of freedom. Discarding those 
particles which bunch, we arrive at 
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where Ni = (^a 2 + ^b 2 + 2ab)~ 1 / 2 , and the correspond- 
ing basis is the same as that in Eq.Q, but the bunching 
terms in |-Bi)li_ri and |-B2)li_ri were discarded. The 
probability of having state Eq.© is (0.5 + ab), it arrives 



at the maximum 3/4 with a = b = 1/2. Eq.Q tells us 
that the bunching is only related to \a p ) and \k p ), and 
these two terms have probability 1/2 of antibunching. 
Thus elements of the density matrix Eq.JSJ differ from 
the initial state upon the normalization factor 1 / \/2 or 
1/2, depending on the total probability of antibunching. 
With the antibunching state as an initial state, next we 
let the second pair of Bob's particles (labelled by 2 in 
the circle) pass through another 50/50 beam splitter and 
perform the same measurement, keeping again the an- 
tibunching results. Having completed the step (4) and 
absorbed all these (n— 1) particles, we get the final state 



= N 2 



( i\) n a 2 






ab 




> |c| 2 




{^Tbc* 


= ((l)"« 2 


f (i) n fe 2 + 2 


= 


\Am\A\) R 


\M ^ 


|4)2|4)£ 


hi) = 


|4)2|4)£ 


1*/) = 


14)214)5 



(72)"- 

\c\ 2 
ab 

( * rbc* 



V2 

2ab)- 1 ' 2 



(I)«c 2 X 

{\Yb 2 ) 



(6) 



and the basis 
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Here, \B t 
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The probability of having this state is Pf - 
{\) n b 2 + 2ab). With n — > oo, the probability pf 



2a6, 



while the entanglement measured by Wootters concur- 
rence is |c| 2 /a6. The entanglement of the final state 

i.e., in the 
2 < ab re- 



reaches its maximal value 1 with |c| = ab 
pure state case. For mixed states, because |c 
quired by Eq.QJ, the entanglement of the final state is 
always less than unity, and it is of relevance to the co- 
herence of state p™ measured by |c| in Eq.Q). On the 
other hand, the state Eq.Q itself is an entangled state, 
the entanglement measured by Wootters concurrence is 
2|c|. To get more entanglement in the final state, it is 
required that |c| > 2ab. Entangled state with |c| < 2ab 
could not be distilled, this is similar to the non-distillable 
entanglement called bound entanglement. 

In conclusion, we have extended the entanglement con- 
centration protocol to the case of mixed states, which 
uses only the effects of particle statistics. The maxi- 
mal entanglement after distillation depends on the off- 
diagonal elements of the initial state that usually is a 
measure of coherence. The distillation efficiency only de- 
pends on the diagonal element of the initial state, it takes 
50% for finite number of particles and tends to 25% for 
infinite large states. 

This work was supported by NCET of M.O.E, and NSF 
of China Project No. 10305002. 



3 



[1] Y. Omar, N. Paunkovic, S. Bose, and V. Vedral, Phys. 

Rev. A 65, 062305(2002). 
[2] N. Paunkovic, Y. Omar, S. Bose, and V. Vedral, Phys. 

Rev. Lett. 88, 187903(2002). 
[3] S. Bose, A. Ekert, Y. Omar, N. Paunkovic, and V. Vedral, 

Phys. Rev. A 68, 052309(2003). 
[4] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 

120(1998). 

[5] A. T. Costa, Jr. and S. Bose, Phys. Rev. Lett. 87, 
277901(2001). 

[6] A. Yu. Kitaev, Ann. Phys. (N.Y.)303, 2(2003). 

[7] M. Plenio and V. Vedral, Contemp. Phys. 39, 431(1998). 

[8] N. Gisin, Phys. Lett. A 210, 151(1996). 

[9] C. H. Bennett, H. J. Bernstein, S. Popescu. and B. Schu- 
macher, Phys. Rev. A 53, 2046(1996). 



[10] C. H. Bennett, D. P. Divincenzo, J. A. Smolin, and W. 

K. Wootters, Phys. Rev. A 54, 3824(1996). 
[11] D. Deutsch etal, Phys. Rev. Lett. 77, 2818(1996). 
[12] S. Bose, V. Vedral, and P. L. Knight, Phys. Rev. A 60, 

1941999. 

[13] L. Hardy and D. D. Song, Phys. Rev. A 62, 052315(2000). 
[14] M. Horodecki, P. Horodccki, and R. Horodecki, Phys. 

Rev. Lett. 78, 574(1997). 
[15] R. T. Thew and W. J. Munro, Phys. Rev. A 63, 

030302(2001). 

[16] P. G. Kwiat, S. Barraza-Lopez, A. Stefanov, and N. 

Gisin, Nature(London) 409, 1014(2001). 
[17] J. W. Pan, C. Simon, C. Brukner, and A. Zeilinger, Na- 

ture(London) 410, 1067(2001). 



